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system at energies for which the underlying classical motion is chaotic. The propagated wavepacket 
is Fourier transformed to yield a scarred eigenstate. 

PACS numbers: 03.65.Sq, 03.65.Ge, 05.45.+b 



Semiclassical mechanics attempts to elucidate and uti- 
lize the relationship between classical dynamics and its 
quantum counterpart in the small Ti regime [1]. For 
bound conservative Hamiltonian systems, two goals may 
"be identified: to develop useful methods for the semi- 
Tnassical propagation of wavepackets and to obtain es- 
fenates for the quantum eigenvalues and eigenfunctions 
Sng classical trajectories. Both have been achieved for 
^d^ssically integrable systems [2]. For classically chaotic 
\^^tems, however, the situation is still unresolved. 
r-HFor the case of eigenfunction or eigenvalue determi- 
nation, most semiclassical studies of bound Hamiltonian 
^^tems have been done in the energy domain. After all, 
rfcirt eigenvalues and eigenfunctions are time independent 
^^antities and should therefore be related to classical 
'ma nifolds that are time invariant. When the correspond- 
iog classical system is integrable, the classical motion is 
^(^tricted to tori and their quantization produces good 
liclassical estimates for the eigenvalues and the eigen- 
ictions [2]. This is in contrast to the case where the 
'T^ssical system is chaotic: the only time invariant struc- 
tires other than the energy shell are then the periodic 
oBjits. Consequently, periodic orbits are at the center of 
fe^ semiclassical efforts in the energy domain for chaotic 
^jifetems [3]. For example, Gutz wilier 's trace formula [4] 
uses the periodic orbits to explain oscillations in the den- 
sity of states and a similar expression due to Bogomolny 
[5] explains the accentuations and attenuations in the av- 
erage coordinate probability density about the classical 
periodic orbits, the so-called scars. 

At present, however, we cannot obtain individual 
eigenvalues or eigenfunctions for generic chaotic Hamil- 
tonians using periodic orbit theory. There are two main 
reasons for this. The first is numerical; to obtain an en- 
ergy resolution off requires all periodic orbits with period 
''' < Tmax ~ 27r7}/e, a quantity which grows exponentially 
in the chaotic regime [6]. Because of this exponential pro- 
liferation of contributing orbits, single eigenvalues can be 
obtained only for very special systems. The second rea- 
son is theoretical and deals with the possible divergence 
of the semiclassical expressions as we try to resolve an 
eigenvalue or eigenstate. This divergence can be traced 
back to the non-commutativity of the limits Ti ^ Q and 



^max ^ CO- 

Similarly, considerable difficulties abound in attempt- 
ing semiclassical propagation in chaotic systems over 
reasonable time periods. For example. Van Vleck- 
Gutzwiller type propagators [2,7] require trajectories sat- 
isfying two-point boundary value conditions. Chaotic 
dynamics makes finding such trajectories extremely dif- 
ficult. Hence, generally useful long time propagators are 
sorely needed to study chaotic dynamics. The avail- 
ability of such methods would also aid in finding de- 
sired eigenvalues and eigenfunctions. In particular, if the 
propagation time is comparable to the Heisenberg time, 
th = 2'iTh/D, where D is the mean level spacing, then 
Fourier transforming the propagated wavepacket exposes 
the semiclassical eigenvalues and eigenfunctions. 

Despite the inherent difficulties, some progress has 
recently been made toward developing useful propa- 
gation methods in the chaotic regime. For example, 
Tomsovic and Heller propagated wavepackets [8], using 
the Van Vleck-Gutzwiller propagator and extracted a 
scarred eigenstate of the stadium billiard [9]. In do- 
ing so they propagated past the so called log time [10], 
iiog ~ j^og Ntf{E) where A is the largest Lyapunov ex- 
ponent of the chaotic fiow [11] and Ntf(E) is the num- 
ber of Planck cells in the phase space volume enclosed 
by a shell of energy E. This time was expected to limit 
the utility of nonuniform propagators of the Van Vleck- 
Gutzwiller type [7] since it was thought to be accompa- 
nied by the exponential proliferation of caustics. How- 
ever, Schulman [12] has recently shown that their success 
in "breaking the log time barrier" arises from patholo- 
gies of the particular system chosen, the stadium bil- 
liard, where most caustics are found near the walls. For 
more realistic systems it is unclear whether non-uniform 
semiclassical expressions can be used for sufficiently long 
times. 

In this letter we use a uniform semiclassical propaga- 
tor to get the time evolution of a wavepacket in a smooth 
Hamiltonian system at energies for which the underlying 
classical motion is chaotic. By uniform we mean that 
there are no caustics (singularities) in our expressions. 
We then Fourier transform this propagation to obtain a 
scarred eigenstate. The method is straightforward and 
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FIG. 1. Contours of the coordinate space probability den- 
sity of the 42nd even (in x) eigenstate, showing striking scar- 
ring by a symmetric libration. The equipotential contour at 
E = 0.4249, the energy of this eigenstate, is also shown. 



FIG. 2. The coordinate representation of the real part of 
the coherent state to be propagated. Solid contours are posi- 
tive, dashed are negative. The contour spacing is 0.2 in Figs. 
2, 3, 4, and 5. 



has no caustics-based time limitations. Furthermore, 
trajectories satisfying two-point boundary values condi- 
tions are not required. 

To obtain the time evolution of a wavepacket |^) we 
consider the following expression: 



(*|e- 



iHt/n 



qi 



rfp/(*|p/)(p/|e-'**/''|q, 



(1) 



Then *(qi,t) = (*|e-'-*(-*)/'*|qi)*. A semiclassi- 
cal approximation for Eq. (1) is obtained by replacing 
{pf\e~*^*'^\qi) with a Van Vleck-Gutzwiller type of ex- 
pression [13]: 



(P/|e-'^*/''|q,).c 



rE- 



(27r;j)"/2^|yd<D)| 

X exp{i(f)/h — iV7r/2} (2) 

The sum is over all classical trajectories that connect q,- 
with pf in a time t; the sum is indexed by the initial 
momenta p,- such that 



q(0) = qi, p(0) = Pi, and p(t) = pf. 



(3) 



The matrix D is one of the four stability matrices which 
describe how trajectories in the immediate vicinity of a 
given classical trajectory behave: 



6q(t) 
6p(t) 



A B\ f6q(0)\ 
C Dj \6p(0)J 



(4) 



The function (t)((li,Pi,t) is a generator of the classical 
motion and is given by 

4>i(ii,Pi,t) = -pitHt) + [ (p(^)q(^) - H)dT. (5) 

Jo 
Finally the index J^[qi,Pi,t] is given by [13,14]: 



'^[qiiPi,^] = — lim arg<^ det(_D -|- ieB) > 

IT e^0+ I- J 



(6) 



with i>[qi,pi,t = 0] = 0. The matrices D and B are de- 
fined by Eq. (4). It can be shown that Eq. (6) is equiva- 
lent to directly derived Maslov indices [15]. 

Inserting the semiclassical expression given in Eq. (2) 
into Eq. (1), we obtain 



(*|e-'**/''|q, 



(27r;i)"/2 



rfp,(*|p(t))|det(i?)|i/2 



X exp{i(f)/h — iV7r/2} (7) 



where p(t) is the momentum resulting from propagation 
beginning at (q_i,Pi). To obtain Eq. (7) we made the 
change of integration variables from p/(= p(t)) to p,-. 
Note that this expression has no singularities. In fact, al- 
though expression (2) develops caustics (i.e. singularities 
caused by the focussing of classical trajectories) as time 
progresses, Eq. (7) is uniform [13]. This expression is in 
the form of an initial value representation [14-16], i.e. it 
allows propagation of trajectories starting at (q,-, p,) and 
does not require [as does, e.g. Eq. (2)] trajectories that 
satisfy two-point boundary conditions. 

We use Eq. (7) to propagate a coherent state (Gaus- 
sian wavepacket) in the following smooth Hamiltonian 
system: 



H = -{pl+pl) + Q-Q^^'+{y-7>^'y 



(8) 



In contrast to billiard systems, this Hamiltonian has a 
smooth potential and is generic in the sense that the 
energy cannot be scaled away. Its classical dynamics have 
been studied extensively [17] and it has been used, in 
the energy domain, to semiclassically study the scarring 
of wavefunctions [18]. The exact quantal wavefunctions 
and energies were obtained as in [18], with k = 0.05. In 
Fig. 1 we show the coordinate probability density of the 
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FIG. 3. The real part of the wavepacket propagated for a 
time t = r/4, where r is the period of the syninietric libra- 
tion. Top: exact quanta! propagation. Bottom: semiclassical 
propagation. 
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FIG. 4. Same as Fig. 3 but for t = r. 



FIG. 5. Same as Fig. 3 but for t = 2r. 

42nd even (in x) eigenstate; it has energy E = 0.4249. 
We also show the simple symmetric libration that scars 

it extensively. 

Here we consider propagation of a wavepacket in this 
classically chaotic regime, extracting this eigenstate semi- 
classically. To do so we we center a coherent state on this 
periodic orbit; the coherent state we chose has the follow- 
ing coordinate representation: 



(q|*) 



1 



1/2 



exp 



(q-q)^ jp 
' 2&2 '^ h 



(q-q/2) 



(9) 



with h = 0.5 and (q, p) is the point where the symmet- 
ric libration intersects the y axis. In Fig. 2 we show the 
real part of the coordinate representation of this coherent 
state, chosen large in coordinate space so that it is well 
localized in momentum space. This feature allows us to 
search out the relevant regions in momentum space where 
Eq. (7) must be evaluated and to carry out the integra- 
tion with low order quadrature methods. Note that we 
do not evaluate Eq. (7) by stationary phase and hence 
avoid the root search problem and caustics associated 
with coalescing trajectories. 

Our calculations of the semiclassical wavepacket prop- 
agation is illustrated in Figs. 3-5 and are compared with 
the exact propagation. The agreement between the exact 
quantal and semiclassical pictures is excellent, even past 



the log time which we estimate to be tiog 



1.5r, where 
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FIG. 6. The Fourier transform of the wavepacket propa- 
gation with E = 0.4249 and tmax = 2.5r. Top: from exact 
quanta! propagation. Bottom: from semiclassical propaga- 
tion. 



r is the period of the symmetric libration. In Fig. 6 we 
show the Fourier transform of our exact and semiclassi- 
cal wavepacket propagations, with Tmax = 2.5r. Since 
we also have the wavepacket propagations backwards in 
time, we expect to get an energy resolution equivalent 
to a propagation time of 5. Or. This is smaller than the 
Heisenberg time, which is th ~ 8. Or at £" = 0.4249. 
However due to special wavepacket we chose, we are able 
to resolve most of the eigenstate and agreement between 
the semiclassical and quantum results is excellent. It is 
important to realize that a similar calculation in the en- 
ergy domain using Bogomolny's formula [18] would not 
produce the nodal structure which is seen to exist along 
the libration, even if all the periodic orbits with periods 
up to Th were included. 

The resultant figures clearly demonstrate that this uni- 
form propagation method is capable of accurate propa- 
gation in chaotic systems over sufficiently long times to 
extract accurate eigenstates. Further studies over longer 
time periods and for larger systems are well warranted 
and are in progress. 



ACKNOWLEDGMENTS 

This work was supported in part by funds provided by 
the Ontario Laser and Lightwave Research Center and 
the U. S. Office of Naval Research under contract number 
N00014-90-J-1014. 



[1] By small ft regime we mean where the dimensionless 
quantity formed from a combination of physical parame- 
ters with h in the numerator can be made small. 
[2] See for example A. M. Ozorio de Almeda, Hamiltonian 
Systems: Chaos and Quantization, Cambridge Univer- 
sity Press (1988); M. C. Gutzwiller, Chaos in Classical 
and Quantum Mechanics, Springer- Verlag (1991); M. S. 
Child, Semiclassical Mechanics with Molecular Applica- 
tions, Oxford University Press (1991). 
[3] See, for example, the following special issue on periodic 

orbit theories: CHAOS 2, (1992). 
[4] M. C. Gutzwiller, J. Math. Phys. 12, 343 (1971); R. 

Balian and C. Bloch, Ann. Phys. (N.Y.) 69, 76 (1972). 
[5] E. B. Bogomolny, Physica D31, 169 (1988); M. V. Berry, 

Proc. R. Soc. London A423 , 219 (1989). 
[6] R. Bowen, Lecture Notes in Mathematics, 470 (1975). 
[7] J. H. Van Vleck, Proc. Natl. Acad. Sci. 14, 178 (1928). 
[8] S. Tomsovic and E. J. Heller, Phys. Rev. Lett. 67, 664 
(1991); M. A. Sepulvida, S. Tomsovic and E. J. Heller, 
Phys. Rev. Lett. 69, 402 (1992). 
[9] S. Tomsovic and E. J. Heller, Phys. Rev. Lett. 70, 1405 
(1993). 

[10] M. V. Berry and N. L. Balasz, J. Phys. A 12, 625 (1979); 
M. V. Berry, N. L. Balasz, M. Tabor, and A. Voros, Ann. 
Phys. (N.Y.) 122, 26 (1979). 

[11] S. Tomsovic, (private communication). 

[12] L. S. Schulman, J. Phys. A 27, 1703 (1994). 

[13] J. R. Klauder, Phys. Rev. Lett. 56, 897 (1986); J. R. 
Klauder, Ann. Phys. (N. Y.) 180, 108 (1987); J. R. 
Klauder, Phys. Rev. Lett. 59, 748 (1987). 

[14] K. G. Kay, Integral Expressions for the Time-Dependent 
Propagator, J. Chem. Phys., in press ; K. G. Kay, Semi- 
classical Propagation for Multidimensional Systems via 
an Initial Value Method, J. Chem. Phys., in press. 

[15] G. Campolieti and P. Brumer, Semiclassical Propagation: 
Phase Indices and the Initial Value Representation, Phys. 
Rev. A., in press. 

[16] W. H. Miller, J. Chem. Phys. 53, 3578 (1970); W. H. 
Miller, J. Chem. Phys. 95, 9428 (1991); G. Campolieti 
and P. Brumer, J. Chem. Phys. 96, 5969 (1992). 

[17] M. Baranger and K. T. R. Davies, Ann. Phys. (N.Y.) 
177, 330 (1987). 

[18] D. Provost and M. Baranger, Phys. Rev. Lett. 71, 662 
(1993); D. Provost, in Coherent States: Past, Present 
and Euture, edited by D. H. Feng, J. R. Klauder, and M. 
R. Strayer, (World Scientific, in press). 



